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Abstract 

We have characterized a novel autosomal recessive Crouzon-like craniosynosto- 
sis syndrome in a 12-affected member family from Antakya, Turkey, the pre- 
senting features of which include: multiple suture synostosis, midface 
hypoplasia, variable degree of exophthalmos, relative prognathism, a beaked 
nose, and conductive hearing loss. Homozygosity mapping followed by targeted 
next-generation sequencing identified a c.479+6T>G mutation in the interleukin 
11 receptor alpha gene (IL11RA) on chromosome 9p21. This donor splice-site 
mutation leads to a high percentage of aberrant IL11RA mRNA transcripts in 
an affected individual and altered mRNA splicing determined by in vitro exon 
trapping. An extended IL11RA mutation screen was performed in a cohort of 
79 patients with an initial clinical diagnosis of Crouzon syndrome, pansynosto- 
sis, or unclassified syndromic craniosynostosis. We identified mutations segre- 
gating with the disease in five families: a German patient of Turkish origin and 
a Turkish family with three affected sibs all of whom were homozygous for the 
previously identified IL11RA c.479+6T>G mutation; a family with pansynosto- 
sis with compound heterozygous missense mutations, p.Pro200Thr and 
p.Arg237Pro; and two further Turkish families with Crouzon-like syndrome 
carrying the homozygous nonsense mutations p.Tyr232* and p.Arg292*. Using 
transient coexpression in HEK293T and COS7 cells, we demonstrated dramati- 
cally reduced ILll-mediated STAT3 phosphorylation for all mutations. Immu- 
nofluorescence analysis of mouse 111 Ira demonstrated specific protein 
expression in cranial mesenchyme which was localized around the coronal 
suture tips and in the lambdoidal suture. In situ hybridization analysis of adult 
zebrafish also detected zfilllra expression in the coronal suture between the 
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Introduction 

The development of the craniofacial skeleton and calvarial 
sutures is a fundamental and complex biological process, 
whose cellular and molecular mechanisms are incom- 
pletely understood. Sutures are key players in the regula- 
tion of the growth and morphogenesis of early skull and 
brain development. Six cranial sutures are present in the 
human cranium: two coronal, one sagittal, two lambdoid, 
and one metopic (Cohen and MacLean 2000). Alterations 
in cellular interactions and signaling between the sutures 
may lead to craniosynostosis, a premature fusion of one 
or more of the cranial sutures during skull development. 
Early suture fusions (single or multiple) lead to distor- 
tions of skull development in the direction of the open 
sutures and are frequently associated with increased intra- 
cranial pressure, impaired blood flow, as well as hearing 
and vision impairments. Moreover, the pressure on the 
growing cerebral cortex results in an increased risk of 
intellectual disability (Renier et al. 1982). 

Craniosynostosis is one of the most common craniofa- 
cial malformations being present in approximately 1 in 
2000 live births (Cohen 1979; Lajeunie et al. 1995) and 
occurs either in rare syndromic, or more frequently, in 
nonsyndromic forms with isolated synostosis. Syndromic 
forms of craniosynostosis are frequently associated with 
limb malformations, but phenotypic associations can 
include a broad spectrum of features with more than 180 
different syndromes associated with craniosynostosis. 
Although the majority of cases are sporadic, approxi- 
mately 10% show familial recurrence (Cohen 2002). 
Monogenic forms of craniosynostosis are usually inherited 



in an autosomal dominant manner with highly penetrant 
mutations in FGFR1-3 and TWIST1 causing the most rec- 
ognizable syndromes including Apert (MIM 101200), 
Crouzon (MIM 123500), Pfeiffer (MIM 101600), Antley- 
Bixler (MIM 207410), Muenke (MIM 602849), and Seat- 
hre-Chotzen (MIM 101400) syndromes. Crouzon syn- 
drome (CS) is characterized by frequent bicoronal 
synostosis and occasional pansynostosis, hypertelorism, 
exophthalmos, divergent strabismus, a beaked nose, short 
philtrum, hypoplastic maxilla, and relative prognathism. 
Malformations of the extremities are more subtle in 
patients with CS than in Pfeiffer and Apert syndromes 
and thus show clinical utility in distinguishing CS from 
other craniofacial syndromes with overlapping cranial 
phenotypes (Kaler et al. 1982; Murdoch-Kinch and Ward 
1997; Mooney and Siegel 2002). Variable inter- and intra- 
familial expressivity of CS is well documented. Typically 
CS is inherited in an autosomal dominant fashion due to 
heterozygous activating mutations in FGFR2 (MIM 
176943). Although autosomal recessive inheritance of CS 
has previously been reported this has received limited 
attention due to the rarity of large families with this 
mode of inheritance (Cross and Opitz 1969; Juberg and 
Chambers 1973). 

Here, we report consanguineous families with a 
Crouzon-like phenotype presenting with multiple suture 
synostosis, exophthalmos, midfacial hypoplasia, and prog- 
nathism without limb malformations. Clinical findings are 
indistinguishable from autosomal dominant CS, although 
intra- and interfamilial variation does exist. Homozygosity 
mapping and targeted next-generation sequencing identi- 
fied missense and nonsense mutations in the IL11RA gene 
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on chromosome 9p21.1-pl3.2 impairing STAT3-related 
downstream signaling. Moreover, we demonstrate that 
mutations in IL11RA also underlie early suture closures in 
pansynostosis. Our data provide exciting evidence for the 
involvement of interleukin 11 signaling in cranial suture 
development and disease. 

Materials and Methods 
Clinical studies 

Three siblings from the index family (Figs. 1A-D and 2B, 
individuals IV:5, IV:6, IV:7) were identified from the 
Hacettepe University Craniomaxillofacial Study Group 
registry. All three affected individuals originated from 
Antakya, Hatay, Turkey, a region with an increased rate 
of consanguinity. A field study was conducted by NAA, 
IV, and SK to evaluate relatives and pedigree construc- 



tion. The complete pedigree structure contained over 427 
individuals and various malformations such as X-linked 
nystagmus (Kaplan et al. 2008), autosomal dominant 
hypodontia, autosomal recessive Carnevale syndrome, and 
craniosynostosis. Only the craniosynostosis cases and their 
families were included in this study. Cases with Crouzon- 
like syndrome were scattered over the various branches of 
this isolate. Affected members, their parents, and survi- 
ving grandparents were examined (Figs. 1 and 2). Blood 
samples were collected and DNA was extracted following 
standard protocols after informed consent was received. 
Institutional ethical board approvals for the research pro- 
ject were obtained. 

Molecular studies 

DNA samples from four affected individuals, their par- 
ents, and unaffected siblings were genotyped using single 




Figure 1. Craniosynostosis phenotypes linked to IL11RA mutations. (A-H) Facial views of representative cases of the Turkish CRS1 family with 
Crouzon-like craniosynostosis. (A and B) Subject IV:5. Facial appearance at 17 years. (C and D) Subject IV: 1 1 at 16 years old. (E and F) Affected 
individual of the CRS3 family carrying the p.Arg292* mutation at approximately 9 years of age. All three subjects were treated by craniotomy at 
an early age. (G and H) Affected individual of the CRS4 family with p.Tyr232* mutation at 12 years of age. (I-K) Skull X-ray with increased digital 
markings and maxillary hypoplasia present in CRS3. 



) 2013 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals, Inc. 



225 



IL11RA and Craniosynostosis 



K. Keupp ef al. 



parents affected individuals healthy individuals 



Chr. 


position (base) 


111:3 111:4 


9 


7036901 





9 


9044621 






9 


11697012 


9 


14030837 




9 


17897828 




9 


20780622 




3 


23721597 


9 


27505967 






32010862 






JHJ1/U44 




9 


72862979 




3 


78021235 




3 


82034241 


Illllll] 

i mil ii 


9 


86662829 


9 


91877285 





IV:5IV:7 IV:1 IV:9 IV:2IV:I0 111:1 I 




D9S319 


29.45 Mb 


D9S1118 


31.82 Mb 


D9S1788 


33.03 Mb 


D9S1845 


33.05 Mb — 




D9S1805 


34.09 Mb 




D9S163 


35.02 Mb 


3. 


D9S1804 


35.86 Mb 




D9S1794 


36.11 Mb 




D9S50 


36.70 Mb — 




D9S1874 


37.12 Mb 



3.4 Mb 



CRS1 family 




-o 



=o 



o 



IV:5 IV:6 



W:7 IV:8 



7> 



=o 



The order of DMA markers 



=o 



D9S1794 

D9SS0 

D9S1874 

1] 



29.45 Mb 
31 .82 Mb 
33.03 Mb 
33.05 Mb 
34.09 Mb 
35.02 Mb 
35.66 Mb 

36.11 Mb 
36.70 Mb 

37.12 Mb 



The order of DNA markers 

D9S319 29.45 Mb 

D9S1118 31.82 Mb 

D9S17B8 33.03 Mb 

D9S1845 33.05 Mb 

D9S1605 34.09 Mb 

D9S163 35.02 Mb 

D9S1604 35.86 Mb 

D9S1794 36.11 Mb 

D9S50 36.70 Mb 

D9S1874 37.12 Mb 



Figure 2. Mapping data of CRS1 family manifesting autosomal recessive Crouzon-like craniosynostosis. (A) Schematic representation of 
homozygosity data of the chromosome 9p21-p12 region. Homozygous genotypes identical to the genotype data obtained from index case IV:5 
(also see pedigree in (B)) are shown in blue. Contrasting homozygous genotypes are shown in white whereas heterozygous SNPs appear orange. 
Noninformativeness as a result of heterozygous genotypes in parent-child trios is indicated in yellow. A single homozygous segment of 
approximately 11 Mb in size between position 27-38 Mb was observed and is marked in the rectangle. The genotyping results of microsatellite 
DNA markers shown on the right are presented in the pedigree below. (B) Pedigree and haplotype analysis of the autosomal recessive CS family. 
Genotyping data and haplotype bars for chromosome 9p markers are shown below the symbol for each individual. Black bars denote the disease- 
associated region. Genotypes in the bracket are the most likely genotypes that are deduced from their children. Thin bars represent 
noninformative genotypes. The critical recombination events positioned the disease allele between DNA markers D9S1845 and D9S50 (marked as 
bold) within a 3.4 Mb critical interval on chromosome 9p21-p12. 
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nucleotide polymorphisms (SNPs) with the GeneChip 
Mapping 10K Array Set (Affymetrix, Santa Clara, CA). 
Genomic DNA (250 ng) was digested by Xbal, followed 
by adaptor ligation and PCR amplification with primers 
provided by the manufacturer (Affymetrix). PCR amplifi- 
cation products were then purified using the Qiagen 
MinElute 96 protocol (Qiagen Inc, Valencia, CA), 
fragmented by DNase I, labeled with terminal deoxynucle- 
otidyltransferase and hybridized to the Mapping 10K Xba 
Chips (Affymetrix). For homozygosity mapping, genotype 
files (CHP files) were generated by the Affymetrix GTYPE 
software and transferred to the VIGENOS (Visual Gen- 
ome Studio) program (Hemosoft Inc, Ankara, Turkey) 
which allows the analysis of genome-wide data in com- 
prehensible visual screens (Kayserili et al. 2009). Haplo- 
types indicating homozygosity by descent were compared 
with those of the index patient IV-5 (Fig. 2A). For 
microsatellite genotyping, several polymorphic DNA 
markers were selected from the critical region on chromo- 
some 9p21-pl3. PCR amplification, followed by denatur- 
ing polyacrylamide gel electrophoresis (6%-7%) and 
silver staining methods were used to separate and visual- 
ize the alleles. Gels were manually photographed (APC 
film; Promega, Madison, WI) and genotyped. The 
MLINK component of the LINKAGE program (FAST- 
LINK, version 3) was used to analyze the linkage data 
(Lathrop and Lalouel 1984; Lathrop et al. 1984; Cotting- 
ham et al. 1993; Schaffer et al. 1994) under the assump- 
tion of an autosomal recessive model with complete 
penetrance. 

Targeted next-generation sequencing 

The entire 3.6 Mb region on chromosome 9 (chr9:33, 
146,801-36,792,821), between markers D9S1845 and 
D9S50, was targeted by array-based sequence capture fol- 
lowed by next-generation sequencing. After stringent 
probe selection by NimbleGen (Roche NimbleGen, Madi- 
son, WI) with uniqueness testing by the Sequence Search 
and Alignment by Hashing Algorithm (SSAHA), a total 
of >2.4 Mb of nonrepeat masked sequences were repre- 
sented on the array, representing 385,000 oligonucleotide 
probes targeting the regions of interest. Sequence capture 
was performed in accordance with the manufacturer's 
instructions (Roche NimbleGen), with the use of the Tita- 
nium optimized protocol as described previously (Niko- 
poulos et al. 2010; Vermeer et al. 2010). In brief, 5 ng of 
proband genomic DNA was used for library preparation 
prior to sequence-capture hybridization. A final amount 
of 3 /ig prehybridization ligation-mediated PCR-amplified 
DNA was hybridized to the customized array, eluted after 
72 h of hybridization and amplified by post hybridization 
LM-PCR. The amplified captured sample was then used 



as the input for emulsion PCR amplification and subse- 
quent sequencing by the Roche 454 GS FLX sequencer 
with Titanium series reagents. 

Mutation screening 

Sequencing of the IL11RA, IL11, IL6, and GP130 genes 
was performed using ABI 310 or ABI 3100 Genetic Ana- 
lyzers (Applied Biosystems, Foster City, CA) following 
standard protocols. All of the identified mutations were 
resequenced in independent experiments, tested for coseg- 
regation with the phenotype within the families, and then 
screened in 100 healthy control individuals from Turkey 
by PCR or restriction digestion (BstXI for the c.710G>C 
mutation, Rsal for c.874 C>T, Hpal for c.696 C>A), 
direct sequencing (c.479+6T>G) or mutation-specific pri- 
mer amplification (c.598C>A). A further 74 control indi- 
viduals from the same village of the original CRS1 family 
near Antakya were screened for c.479+6T>G by direct 
sequencing analysis. Primer sequences and conditions are 
provided in Table SI. 

Patient cDNA and in vitro exon-trapping 
analyses 

Patient RNA was extracted from fresh blood using the Pax- 
gene Blood RNA system (Qiagen, Hilden, Germany). RT- 
PCR was performed via RevertAid First Strand cDNA syn- 
thesis Kit (Fermentas, St. Leon-Rot, Germany). RT-PCR 
products were amplified by standard PCR using IL11RA 
specific primers, visualized by agarose gel electrophoresis 
and used for direct sequencing. Primers for amplification 
and sequencing were designed according to the reference 
sequences for exons 2 and 7 of IL11RA. In vitro analysis of 
splice mutation c.479+6T>C was performed by the Exon- 
Trapping System (Gibco Invitrogen, Karlsruhe, Germany). 
Genomic fragments of the wildtype and the mutant variant 
of IL11RA including exon 4 and 5 and the flanking intron- 
ic sequences were cloned into the pSPL3 expression vector 
using the EcoRI and EcoRV restriction sites. HEK293T 
cells were transiently transfected with wild type and 
mutant constructs followed by RNA extraction using 
RNeasy Kit (Qiagen) and reverse transcription via Revert- 
Aid First Strand cDNA synthesis Kit (Fermentas). 

IL11RA expression constructs 

The pSVL-sILllRflag vector described previously (Pflanz 
et al. 1999), was used as the template for the generation 
of IL11RA expression constructs. The human cDNA of 
the soluble isoform of IL11RA fused to a C-terminal Flag 
tag was amplified via PCR and cloned into the pRK5 
expression vector using NotI and Hindlll restriction sites. 
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The identified mutations (p.Pro200Thr, p.Arg237Pro, 
p.Tyr232*, and p.Arg292*) were introduced by site-direc- 
ted PCR mutagenesis with primers containing the specific 
nucleotide substitutions. 

Cell culture and transient transfection 

HEK293T and COS7 cells were cultured in Dulbecco's 
Modified Eagle Media (DMEM) supplemented with 10% 
fetal bovine serum (FBS; Gibco) and antibiotics. Cells 
were transiently transfected with vectors containing wild 
type and mutant variants of IL11RA cDNA or genomic 
fragments of IL11RA using Lipofectamine 2000 
(Invitrogen) following the manufacturer's instructions. 

Western blot analysis 

Eighteen hours after transfection, cells were activated by 
the addition of 0.5 fig IL11 (Peprotec, Rocky Hill, NJ) 
into 8 mL growth medium for 10 minutes at 37°C. After 
incubation, activated and nonactivated control cells were 
harvested and total proteins solubilized using ice-cold 
lysis buffer. The total protein concentration of extracts 
was determined using the BCA Protein Assay Kit (Pierce 
Protein Research Products; Thermo Fischer Scientific, 
Rockford, IL), and 15 jig of total protein from each sam- 
ple was separated by 4-12% SDS-PAGE and blotted onto 
nitrocellulose membranes (GE Healthcare, Mtinchen, Ger- 
many). Protein detection was performed using antibodies 
against pStat3 (Santa Cruz Biotechnology Inc, Santa 
Cruz, CA). Equal protein loading and transfection effi- 
ciency was confirmed by reprobing of the membranes 
with antibodies against Flag and /?-Actin. Peroxidase con- 
jugated secondary antibodies were purchased from Santa 
Cruz and blots were developed using an enhanced 
chemiluminescence system, ECL Plus (Amersham, UK), 
followed by exposure on autoradiographic film (GE 
Healthcare). 

Expression analysis 

For histological analysis, P0 wild-type mouse heads were 
embedded in tissue freezing medium (Leica Biosystems, 
Nussloch, Germany) and frozen by immersion in isopen- 
tane precooled in liquid nitrogen. Twelve micrometer cryo- 
sections were stained with hematoxylin and eosin (Merck, 
Darmstadt, Germany). For immunofluorescence analysis, 
12 f,im cryosections were blocked in 1 x PBS containing 
0.5% TritonX and 10% Fetal Calf Serum and probed with 
an antibody raised against 111 Ira (1:20, (N-20) sc-993, 
Santa Cruz). Antibody binding was visualized with a 
Cy3-conjugated secondary antibody (Invitrogen) and the 
sections counterstained with DAPI (Invitrogen). 



In situ hybridization 

In situ hybridizations were performed using probes gener- 
ated from cloned cDNA fragments of zebrafish ill Ira. 
Probes were synthesized from linearized plasmids using 
the Roche digoxigenin RNA synthesis kit. Adult zebrafish 
(6-12 months old) were fixed in 4% paraformaldehyde, 
decalcified in 0.5 mol/L EDTA (pH 7.4) for 5 days, dehy- 
drated, cleared, embedded in paraffin, and sectioned at 
8 jim. In situ hybridizations were performed as previously 
described (Moorman et al. 2001) with the following 
modifications: prehybridization and hybridization were 
performed at 62°C and BMPurple (Roche, Mannheim, 
Germany) was used for chromogenic detection. 

Results 

The CRANIRARE consortium (http://www.cranirare.eu) is 
an E-RARE funded network and represents an integrated 
clinical and molecular approach for craniofacial malfor- 
mations aiming to identify the molecular pathogenesis for 
various craniofacial malformations. The strategy to iden- 
tify novel genes within the presented study was based pri- 
marily upon a large collection of patients from the 
Hacettepe University Craniomaxillofacial Study Group 
registry, which principally includes syndromic forms of 
craniosynostosis. A prescreen for mutations in the crani- 
osynostosis genes FGFR1, FGFR2, FGFR3, and TWIST was 
performed and mutations identified in approximately 
50% of cases (B. Wollnik, pers. comm., 2010). 

An autosomal recessive Crouzon-like 
phenotype 

The index family (Fig. 1) was identified due to CS-like 
findings segregating with an autosomal recessive mode of 
inheritance. Consistent findings were multiple suture syn- 
ostosis, midfacial hypoplasia, variable exophthalmos, and 
relative prognathism (Table 1). A beaked nose and con- 
ductive hearing loss are features observed in some family 
members. The index case (Figs. 1A and B, 2B individual 
IV:5) is a 17-year-old female with complex suture synos- 
tosis, brachycephaly, and hearing loss with an initial diag- 
nosis of CS. Her twin sisters (Fig. 2B, individuals IV:6 
and IV:7) had similar craniofacial malformations while 
both parents were unaffected. The index case was treated 
with a strip craniectomy at the age of 3 years and her 
twin sisters had cranial surgery at the age one. A total of 
twelve individuals with Crouzon-like phenotype were 
identified in various branches of the pedigree (Fig. 3A 
and Table 1) following practitioner review. Two individu- 
als (Fig. 3A, IV: 11 and 111:8) had dental extractions due 
to failure of tooth eruption. 
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Figure 3. Identification of a donor splice-site mutation in IL11RA. (A) Pedigree of the large Turkish kindred with Crouzon-like craniosynostosis. 
Capital letters beneath each symbol describe the genotypes of family members for the two variations identified by targeted next-generation 
sequencing in IL11RA and PIGO. Amino acid and nucleotide substitutions are shown in red. (B) Sequence chromatograms showing the identified 
c.479+6T>G mutation in IL11RA in the affected individual IV: 1, heterozygous carrier 111:1, and a wild-type (wt) control. The arrow indicates the 
nucleotide substitution on the sixth position behind the exon intron boundary. (C) mRNA analysis of IL11RA transcripts in an affected individual 
carrying the c.479+6T>G mutation in the homozygous state. RNA was reverse transcribed and after PCR amplification separated on a 1 % 
agarose gel resulting in two visible bands containing three fragments (patient lane). The upper band (-600 bp) contained a WT fragment and an 
aberrant transcript which includes additional 5 bp generated by the use of an alternative donor site introduced by the mutation. The lower band 
(-450 bp) represents an aberrant transcript generated by skipping of exon 4. Frequencies of transcripts are indicated at the right and represent 
the analysis of 21 clones after subcloning and sequencing of fragments. 



Identification of a homozygous splice-site 
mutation in IL11RA 

A total of 12 individuals from this family were genotyped 
using the Affymetrix GeneChip Human Mapping 10K 
SNP Array. Genotype files were generated with Affymetrix 
GTYPE software and transferred to the VIGENOS (Visual 
Genome Studio) program (Hemosoft) (Kayserili et al. 
2009) to identify the minimum region of shared homozy- 



gosity (Fig. 2A). A single homozygous region on chromo- 
some 9p21-pl2 was found in all the affected individuals, 
for which the parents and healthy siblings were heterozy- 
gous. Subsequently, typed microsatellite markers for this 
chromosomal regions included additional family members 
(Fig. 2B). Homozygosity was confirmed and recombina- 
tions between markers D9S1845 and D9S50 defined a crit- 
ical interval of 3.4 Mb containing 78 annotated genes. 
Linkage analysis obtained a maximum pairwise lod score 
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Table 2. Next generation sequencing variant statistics. 



Filter applied 


Number of variants 


All variants called 


2529 


Mapping to target region on chr.9 


1958 


Not previously known SNP (dbSNP129) 


105 


Homozygous variant (>80% variant reads) 


74 


Coding/splice site variant 


2 



of 2.306 at Theta = 0 cM between the disease allele and 
the DNA marker D9S1805 (data not shown). Subse- 
quently, array-based sequence capture of the complete 
3.4 Mb region was performed, followed by next-genera- 
tion sequencing on a one-quarter plate of a Roche 454 
sequencing run. Of the 61 million sequencing reads 
approximately 50% of the sequence data mapped back to, 
or near, the targeted region of chromosome 9p21-pl2 
and the average of coverage for the region was 11 -fold. 
98.4% of all targeted bases were sequenced at least once 
and >83% of all targeted bases were covered at least 
fivefold. 

Among the identified homozygous variants (Table 2) 
only two were not annotated as SNPs and were located 
within either exonic sequences or splice sites: a nonsynon- 
ymous coding variant, p.Gly338Gln, in PIGO and the 
c.479+6T>G variant of the donor splice site of exon 4 of 
IL11RA. Both variants were confirmed by Sanger sequenc- 
ing, but only the c.479+6T>G variant (Fig. 3B) cosegre- 
gated with the disease in the extended pedigree and was 
present in all 12 affected individuals in various branches 
of the family (Fig. 3A). The c.479+6T>G mutation was 
not present in 100 Turkish controls. Further evidence that 
this splice-site mutation is likely to be the pathogenic 
allele came from the finding that this mutation leads to a 
high percentage of aberrant IL11RA mRNA transcripts in 
an affected individual (Fig. 3C) and alters normal splicing 
in an in vitro exon-trapping experiment (Fig. SI). It is 
noteworthy that the donor splice site of exon 4 is not 
completely abolished and approximately 25% of the nor- 
mal transcript was observed in mRNA derived from one 
affected individual. 

In order to determine the carrier frequency of the 
IL11RA c.479+6T>G mutation, we tested a total of 74 
individuals from the same village. From these data, we 
included parents and spouses only (30 individuals) to 
calculate the mutant allele frequency. Expected carrier 
frequency (2pq) was approximately 0.30 for this village 
(normal T-allele frequency (p) = 0.816 and mutant 
G-allele frequency (q) = 0.184). Following identification 
of the IL11RA c.479+6T>G mutation, we confirmed the 
presence of this allele in an additional consanguineous 
family with three affected siblings with Crouzon-like find- 
ings from a nearby village of Antakya. Although no clear 



relationship between this and the index family could be 
established, the identical c.479+6T>G mutation was iden- 
tified by direct sequencing. In conclusion a total of 15 
affected individuals with the IL11RA c.479+6G>T muta- 
tion suggests a founder effect in the Antakya region of 
Turkey. 

Additional IL11RA mutations in Crouzon-like 
syndrome and pansynostosis 

All 13 coding exons of the IL11RA gene were screened by 
DNA sequencing in a cohort of 79 patients with a diagnosis 
of either CS, pansynostosis, or nonclassified syndromic 
craniosynostosis. The IL11RA c.479+6T>G mutation was 
also identified in a consanguineous patient from Germany 
of Turkish origin. Interestingly, this patient presented with 
pansynostosis suggesting that mutations in IL11RA may 
display clinical variability (Table 1). Pansynostosis was also 
the leading diagnosis in two affected siblings from a non- 
consanguineous family also from Turkey (Table 1), both of 
whom had the following compound heterozygous missense 
mutations (Fig. 4A): IL11RA c.598C>A (p.Pro200Thr) 
paternal and c.710G>C (p.Arg237Pro) maternal. Both 
mutations are located within the FN3 domain of the 
IL11RA protein (Fig. 4B,C). The homozygous nonsense 
mutations, c.696C>A (p.Tyr232*) and c.874C>T 
(p.Arg292*), were identified in two additional Turkish 
families with Crouzon-like syndrome (Table 1 and 
Fig. 4A). These mutations are predicted to truncate the 
IL11RA protein within the FN3 domain, leading to a loss of 
the C-terminal transmembrane domain (Fig. 4B,C). Nei- 
ther of these mutations were found in 100 healthy Turkish 
control individuals. 

Impaired IL11RA signaling caused by 
identified mutations 

IL11RA is a known receptor involved in interleukin sig- 
naling. After the binding of its ligand, IL11, the receptor 
binds to the cytokine signal transducer GP130, leading to 
dimerization and subsequent transphosphorylation of 
intracellular tyrosines within the cytoplasmic domain 
through the associated Janus kinases (Jaks) (Dahmen 
et al. 1998). One of the known downstream targets is 
STAT3, which is activated by phosphorylation in the sig- 
naling cascade (Fig. 5A) (MIM 102582). Using transient 
expression experiments in HEK293T and COS7 cells, we 
demonstrated that both the missense mutations and non- 
sense mutations of IL11RA reduced in vitro IL11- 
dependent STAT3 phosphorylation to undetectable levels 
(Fig. 5B). These results clearly implicate impaired inter- 
leukin 11 signaling as a novel pathogenic mechanism 
underlying Crouzon-like syndrome and pansynostosis. 
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Figure 4. Additional IL11RA mutations. (A) Pedigrees and sequence chromatograms of additional IL11RA mutations. Capital letters above the 
pedigrees describe the amino acid changes found in affected family members. Substitutions are marked in red. (B) Schematic representation of 
the genetic and molecular structure of IL1 1 RA and positions of identified mutations within IL1 1 RA domains. (C) Predicted protein structure of the 
IL11-IL11RA receptor complex. Amino acids at the positions of the identified missense (Pro200, Arg237) and nonsense mutations (Tyr232, 
Arg292) are shown. 



Mutation screen in IL11RA pathway genes 

In order to test the hypothesis that mutations in addi- 
tional genes encoding proteins involved in IL11RA signal- 
ing could also be associated with craniosynostosis, the 
five coding exons of IL11 were sequenced in 46 patients 
with various forms of craniosynostosis as well as the five 
coding exons of IL6 in 28 patients and the 14 coding 
exons of the GP130 in 44 patients. No candidate causative 
mutations were identified, from which we conclude that 
alterations of these genes are not a frequent cause of 
syndromic craniosynostosis. 

Conserved expression of 1111 ra in mouse 
and zebrafish sutures 

Immunofluorescence analysis of mouse 111 Ira demon- 
strated specific protein localization in cranial mesen- 
chyme and in the granular layer of the epidermis of 
newborn mice (Fig. 6A,B). Histological analysis of PO 



mouse cranial tissue revealed the normal structure of the 
coronal suture, formed by the overlapping frontal and 
parietal bones (Fig. 6A, upper panel). Immunofluores- 
cence analysis of the murine ortholog of IL11RA demon- 
strated protein localization around the coronal suture 
tips (Fig. 6A, lower panel). Significant ILllra protein 
levels could also be detected in the lambdoid suture in 
the posterior skull (Fig. 6B). Here, the suture is not 
formed by overlapping plates but by direct contact of 
the two opposing bones. Ill Ira protein could be detected 
in the zone of bone contact between the plates forming 
the lambdoid suture (Fig. 6B). In situ hybridization anal- 
ysis of adult zebrafish also confirmed cranial expression 
of the zebrafish ortholog of IL11RA (Fig. 6C). The zebra- 
fish ill Ira RNA was localized in the coronal suture, 
between the overlapping frontal and parietal plates 
(Fig. 6C, left panel). The arrangement of the skull plates 
in fish is very similar to that observed in mouse (Quarto 
and Longaker 2005). In the anterior part of the zebrafish 
head, two further joint-like structures of flat bones were 
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Figure 5. Impaired IL1 1RA signaling causes craniosynostosis. (A) Simplified schematic view of the signaling pathway involving IL1 1RA. Binding of 
IL1 1 to IL11RA is thought to trigger homodimerization of the GP130 receptor and transphosphorylation of tyrosines of GP130 through the 
associated JAKs. Subsequently, STAT3 is activated by tyrosine phosphorylation. (B) Western blot analysis of IL11 -induced phosphorylation of 
STAT3 in transiently transfected HEK293T and COS7 cells. Cells were activated with IL1 1 (+) or left unactivated (-). Transfection efficiency and 
equal protein amounts were confirmed by reprobing of the membrane with antibodies against Flag and B-Actin. C-terminal fused Flag-tag was 
neither detectable in proteins with premature stop signals nor in untransfected control cells. 



detectable and specific ill Ira expression was also detected 
in the areas of these bone contacts (Fig. 6C, right panel). 
The observed similarities between mouse and zebrafish 
cranial structure provides evidence for a conserved 
mechanism of suture formation in vertebrates. Further- 
more, we can conclude a conserved involvement of 
IL11RA in suture development in both zebrafish and 
mouse. 

Discussion 

Here, we demonstrate that mutations in the IL11RA gene 
cause autosomal recessive Crouzon-like craniosynostosis 
or syndromic pansynostosis with mild craniofacial pheno- 
type via impaired interleukin- 1 1 mediated signaling, based 
on mutations in five families. 

Consistent findings in all affected family members of 
the CRS1 family were multiple suture synostosis, midfa- 
cial hypoplasia, variable exophthalmos, prognathism, a 
variably beaked nose, and conductive hearing loss. The 
craniosynostosis observed was complex and affected most 
sutures (pansynostosis) in these individuals with an initial 
diagnosis of CS. CS is a well-known entity and has been 
associated with heterozygous gain-of-function mutations 
in the FGFR2 receptor (MIM 176943). It is mainly char- 
acterized by frequent bicoronal synostosis with occasional 
pansynostosis, hypertelorism, exophthalmos, divergent 
strabismus, beaked nose, short philtrum, hypoplastic 



maxilla, and prognathism. This study confirms previous 
reports (Cross and Opitz 1969; Juberg and Chambers 
1973) that an autosomal recessive form of a Crouzon-like 
craniosynostosis exists and shows that this can be caused 
by mutations in IL11RA. After initial gene identification, 
we screened all syndromic and nonsyndromic craniosy- 
nostosis cases in the CRANIRARE database and detected 
mutations only in cases classified as having CS. Therefore, 
we suggest that patients with Crouzon-like features and 
parental consanguinity should be first tested for IL11RA 
mutations. No additional skeletal anomalies affecting the 
cervical spine, wrist and limbs as well as no internal organ 
involvement was observed. We did not observe specific 
phenotypic differences between IL11RA and FGFR2 
caused craniosynostosis. 

Nieminen et al. (2011) published a study showing that 
recessive mutations in IL11RA cause craniosynostosis, 
delayed tooth eruption, and supernumerary teeth. The 
authors state that one of the children in Family 2 had a 
Crouzon-like facial appearance which was not considered 
typical of any of the previously described craniosynostosis 
syndromes. In accordance with the clinical data presented 
by Nieminen et al, we observed class I and class II mal- 
occlusion and failure of tooth eruption in a large number 
of patients (Table 1). These findings are also consistent 
with IL11RA having a role in dental development. 

It is known that heterozygous activating mutations in 
FGFR2 can cause different phenotypes, such as Crouzon, 
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Figure 6. 111 Ira expression in mouse and zebrafish sutures. (A-B) Sagittal cryosections of PO mouse heads; scale bars: 100 nm. (A) Upper row, 
histological staining with hematoxylin and eosin shows the coronal suture (arrow heads) and surrounding tissue. Lower row, immunofluorescence 
staining displays 111 1 ra localization (red) in the granular layer of the skin (arrow) and around the coronal suture tips (asterisks). Nuclei were stained 
with DAPI (blue). White lines mark cranial bone plates. (B) Immunofluorescence staining of 111 1 ra (red) shows protein in the lambdoidal suture 
(asterisks). Nuclei were stained with DAPI (blue). White lines mark cranial bone plates. (C) In situ hybridization of H11ra mRNA on sagittal paraffin 
sections of adult zebrafish heads. H1 1ra is expressed within the coronal sutures (left) and anteriorly in craniofacial bone fusions (right, arrows) 
within areas of bone contact; scale bars: 100 /im. 



Apert, Pfeiffer, or Antely-Bixler syndromes. The factors 
modifying the phenotypic outcome of FGFR2 mutations 
are not well understood. Interestingly, we provide evi- 
dence that recessive mutations in IL11RA also show vari- 
able clinical expression (Table 1). As an example, both 
affected individuals of the CRS2 family presented with 
only pansynostosis, that is, without midfacial hypoplasia, 
exophthalmos, or relative mandibular prognathism, indi- 
cating that the clinical spectrum of IL11RA mutations 
ranges from pansynostosis to Crouzon-like craniosy- 
nostosis. The pansynostosis in this family is caused by 
compound heterozygous missense mutations, p.Pro200- 
Thr and p.Arg237Pro, in IL11RA. Functional analysis of 
mutant proteins in two different cell systems clearly indi- 
cates that both mutations cause an impairment of 
IL1 IRA-mediated signaling as evidenced by the lack of 
downstream STAT3 activation (Fig. 5B). No functional 
differences in STAT3 activation were observed between 
the missense and nonsense mutations analyzed. Therefore, 
it is unlikely that variable phenotypes observed in our 



patients can be attributed to the different genotypes 
observed. This is further supported by the fact that the 
affected individual of the CRS5 family, who carries the 
homozygous c.479+6T>G mutation, which is also present 
in 15 Crouzon-like craniosynostosis cases from the region 
of Antakya, does present with Crouzon-like syndrome but 
has pansynostosis. It will be an important aim for future 
studies to identify molecular modifiers responsible for 
phenotypic variability. 

Recent advances in massive parallel sequencing tech- 
nologies have accelerated gene identification studies (Ng 
et al. 2010; Kalay et al. 2011). Here, targeted NGS of 
the complete linked critical region was used to sequence 
all 78 annotated genes within this region. Interestingly, 
only two variations of a putatively causative nature 
were detected (Table 2) of which only one cosegregated 
with the phenotype in the extended pedigree. These 
data nicely underscore the power of NGS applications, 
especially when mapping information in a family is 
available. We identified one homozygous splice-site 
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mutation (c.479+6T>G), two homozygous nonsense 
mutations (p.Tyr232* and p.Arg292*), and compound 
heterozygosity for two missense mutations (p.Pro200Thr 
and p.Arg237Pro) in our patient cohort. All mutations 
identified in our study were shown to cause impairment 
of ILllPvA function upon activation by IL11. We dem- 
onstrated that the c.479+6T>G donor splice-site muta- 
tion causes incomplete aberrant mRNA splicing, and 
that both missense and nonsense mutations fail to lead 
to STAT3 phosphorylation in response to IL11. There- 
fore, we conclude that impaired IL11-IL11RA signaling 
is the pathogenic mechanism underlying autosomal 
recessive Crouzon-like craniosynostosis and autosomal 
recessive pansynostosis. No mutations were identified in 
the IL11, IL6, and GP130 genes in the CRANIRARE 
cohort. 

Interleukin- 1 1 (IL11) belongs to the interleukin-6-fam- 
ily of cytokines (Paul et al. 1990; Yin et al. 1992) and has 
mesenchymal cell specific activity, which includes bone 
marrow stromal cells and osteoblasts (Girasole et al. 
1994). IL11 has a wide spectrum of biological functions 
and is involved in hematopoiesis, immune responses, 
neuronal processes, and bone metabolism (Musashi et al. 
1991; Yin et al. 1992; Mehler et al. 1993; Du and 
Williams 1994, 1997). Ill Ira plays a role in bone forma- 
tion and remodeling and is expressed in both osteoblasts 
and mature osteoclasts. Moreover, bone-forming and 
bone-resorbing cells are potential targets of 111 1 (Romas 
et al. 1996). During bone developmental processes, cyto- 
kines, including IL11, which binds to gpl30 receptor, can 
modulate stromal cell commitment to osteoblast differen- 
tiation pathways (Gimble et al. 1994). IL6 induced osteo- 
blast maturation is mediated primarily by activation of 
the IAK/STAT signaling pathway (Bellido et al. 1997). 
IL11 also stimulates RANKL expression in osteoblasts and 
stromal cells, thereby indirectly promoting RANK-medi- 
ated osteoclastogenesis (Yasuda et al. 1998). Furthermore, 
IL11 signaling has direct effects on osteoclast precursors 
and stimulates osteoclast formation in a RANKL-indepen- 
dent mechanism (Kudo et al. 2003; Sims et al. 2005). 
Both direct and indirect ILll-stimulated osteoclast forma- 
tion is a gpl30-dependent process. However, RANKL- 
induced indirect osteoclast promotion occurs via a 
STAT3-mediated process (O'Brien et al. 1999), whereas 
direct effects on osteoclast formation are mediated by 
SHP2/Ras/MAPK signaling (Sims et al. 2004, 2005). Mice 
deficient for 111 Ira (ttllml~/~) exhibit abnormal bone 
turnover and long bone structure and moreover, a reduc- 
tion of osteoclast differentiation has been shown to be 
responsible for increased trabecular bone volume (Sims 
et al. 2005). Recently, disturbed cranial growth and suture 
activity in III Ira knockout mice have also been described 
(Nieminen et al. 2011). In this study, expression analysis 



in newborn mice indicated a highly specific localization 
of 111 Ira within the mesenchyme of the coronal suture. 
Interestingly, during suture formation in adult zebrafish, 
illlra transcripts were also detectable in the intervening 
mesenchymal tissue, which provides further evidence for 
a highly conserved role of IL11RA during suture 
development in vertebrates. 

The cranial skull vault anatomy and development is 
very similar in mouse and fish. For both mouse and ze- 
brafish it is described that sutures - except the posterior 
frontal suture in mice - remain patent for the life of 
the animal (Opperman 2000; Quarto and Longaker 
2005). The highly specific detection of Illlra restricted 
to both ends of the cranial plates within the suture 
matrix in mice and expression of illlra in a fibrous 
stripe of cellular material within the zebrafish coronal 
sutures reveal an involvement of Illlra in separating the 
frontal and the parietal bones, which subsequently could 
prevent a fusion of these bones. This suggestion is sup- 
ported by these genetic and clinical findings showing 
that mutations in IL11RA cause premature fusion of cra- 
nial sutures in patients with Crouzon-like syndromes 
and pansynostosis. 

In our mouse expression study, Illlra localized to a 
specific cell cluster within the mesenchyme, possibly rep- 
resenting a functionally distinct cell population. To date, 
it is not clear if alterations in Illlra cause a dysfunction 
of osteoblasts or osteoclasts or an early developmental 
effect on mesenchymal precursor cells during osteoblast/ 
osteoclast differentiation. Hyper-ossification defects as 
evident in the fused sutures of the patients are most 
likely due to increased osteoblast or decreased osteoclast 
activity. However, direct STAT3-mediated IL signaling, 
which according to our in vitro data is compromised by 
the IL11RA mutations, has been shown to have a stimu- 
latory effect on osteoblasts, whereas direct stimulatory IL 
signaling on osteoclasts is mediated by other transduc- 
tion pathways (see above). Therefore, we hypothesize 
that the pathogenesis of the disorder in IL11RA patients 
is caused by a disruption of the osteoblast-mediated, 
RANKL-dependent, indirect positive effect on osteoclas- 
togenesis or osteoclast activity. However, other effects 
on osteoblasts independently of their regulatory role on 
osteoclasts also cannot be ruled out, such as an attenu- 
ated recruitment of osteoblasts from mesenchymal pre- 
cursors during earlier development and a corresponding 
aberrantly high number of osteoblasts at the relevant 
stages. 

In conclusion, our results provide evidence for a crucial 
and conserved role of IL11RA during craniofacial devel- 
opment and suture formation. We propose an inhibitory 
effect of Illlra within sutures, thereby preventing their 
premature fusion, something which was also recently sug- 
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gested in Illlral~/~ mice. In this study, it was shown 
that IL11 signaling has a negative regulatory effect on 
suture closure (Nieminen et al. 2011). Taken together, 
our data demonstrate a vital role for IL11-IL11RA signal- 
ing during craniofacial morphogenesis in humans, mouse, 
and zebrafish. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Fig. SI. In vitro exon-trapping analysis in HEK293T cells. 
(A-B) Transiently expressed splice vector pSPL3, which 
contained a genomic fragment of IL11RA including exon 
4 and 5 and surrounding intronic sequences, resulted in 
two alternative splice products. Product 1 (upper band) 
contained a fragment with correctly spliced exon 4 and 5 
of IL11RA, and was only present in cells transfected with 
WT constructs. In product 2 (lower bands) exon 4 is 
skipped and was present in cells expressing WT and 
mutant constructs containing splice mutation 
c.479+6T>G. 

Table SI. Primer sequences used for amplification and 
sequencing of the IL1 IRA gene from genomic DNA. 
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